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Abstract 

We analyzed the atmospheric neutrino experiments of SuperKamiokande 
including zenith angle dependence's using the three-flavor neutrino framework 
with the hierarchy mf ~ tti-I ^ m^. Taking into account the terrestrial, 
solar neutrino experimental data and the atmospheric neutrino experiments 
including the sub-GeV and multi-GeV data in SuperKamiokande, large angle 
solution in the solar neutrino experiments is favored and the range of the 
mass parameter Am23 restricted between O.OSeV^ — 2eV'^. Allowed regions 
of mixing parameters are (6*13 < 4°, 27° < ^23 < 32°) for Am23 = leV^ and 
(6*13 < 3°, 28° < 6*23 < 33°) for Aml^ = O.leV^. 



I. INTRODUCTION 



Observing the neutrino oscillations is one of the most important experiments with respect 
to new physics beyond the Standard Model. In recent SuperKamiokande atmospheric experi- 
ments [0,0 , definite atmospheric neutrino anomaly is observed following that in Kamiokande 
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experiments and zenith angle dependence's of double ratios R{fi/e) in atmospheric neu- 
trinos are obtained, those which have decreasing characters of double ratios to increasing 
of zenith angle for both sub-GeV and multi-GeV experiments. Decreasing character sug- 
gests the strong evidence for neutrino oscillations and observing these decreasing characters 
in both sub-GeV and multi-GeV experiments give definite information about the neutrino 
mixing parameters; Am^ and sin^ 29. 

In the terrestrial neutrino experiments, many experimental groups have tried to observe 
the neutrino oscillations. The E531 ||^, CHORUS and NOMAD experiments using the 
beam from accelerators search for u^- appearance in u^, and E776 KARMEN ^ and 
LSND IP] experiments using the accelerator beams are searching for v^, and z/^ — > 

oscillations. The experiments using nuclear power reactor [^jO] search for the disappear- 
ance of z/g, in which z/g — »■ vx{X = fi, r) transitions are expected. These experiments do not 



observe significantly large neutrino transitions. Especially, the recent Chooz experiment |Tl 
excludes the parameter region given approximately by Am^ > 0.9 x 10~'^eV^ for maximum 
mixing and sin^ 29 > 0.18 for large Am^. 

In this paper, we analyze the atmospheric neutrino experimental data obtained in Su- 
perKamiokande including the zenith angle dependence in three-flavor neutrino framework 
with a hierarchy of neutrino masses mi ^ m2 <^ m^. In previous work |]12[, we have analyzed 
the solar neutrino, terrestrial neutrino and atmospheric neutrino experiments in three-flavor 
neutrino framework with hierarchy mi ^ m2 -C including the zenith angle dependence's 
partially. In this work, we perform the full analysis of zenith angle dependence's in Su- 
perKamiokande atmospheric neutrino experiment |0] and analysis of Chooz data |]11| in 
terrestrial experiment. 

In three-flavor neutrino framework with hierarchy mi ^ m2 <^ ms, there are 5 parameters 
sin^2^i2, 6'i3, 6*23, Am^2 = ""^2 ~ '^"^ia ~ ""^s ~ ""^2 concerned with the neutrino 
oscillation. For solar neutrino deficit, MSW solution [O^U 'Hi considering the matter effects 
predicts the large angle solution sin^ 2^i2 = 0.6—0.9, Ami2 = 4x 10^® — 7x lO^^eV^ and small 
angle solution sin^ 2^12 = 0.003 - 0.01, Amj^ = 3 x lO"*' - 1.2 x lO^^eV^ for ^13 = 0° - 20°, 
and these large angle and small angle solution are merged for ^13 = 25° — 50°. The vacuum 



solution in solar neutrino is obtained as Ami2 ~ 10 ^"^eV^ |T5|. In terrestrial neutrino 



experiments, Fogli et al. ||T6[| have analyzed thoroughly and obtained the allowed regions 
on tan6'J3 — tan^ 6'23 plane for various values of mass parameter Am23. We also analyzed 
the terrestrial neutrino experiments in which we have not included Chooz experiment, 
and obtained the similar results as the one of Fogli et al. Obtained allowed region of 
tan^6'i3 — tan^6'23 plane, for example for Am23 = 0.02eV^, is tan^ 6*13 > 25,tan^6'i3 < 0.04 
for all values of tan^ 6*23. For higher values of Am23, allowed region on tan^ 6*13 — tan^ ^23 
plane was more restricted, and for lower values of it, the allowed region spread out and for 
Am23 < O.OleV^ there was no restriction on tan^ ^13 — tan^ ^23 plane. 

Atmospheric neutrino experiments have been analyzed in three-flavor neutrino framework 
with neutrino mass hierarchy mi ~ m2 -C by many authors |1^|T^JT^ . The atmospheric 
neutrino anomaly of Kamiokande experiments |]TU[ has been explained by ^ oscillation 
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with a parameter Am^^ ~ 10~^ — 10~^eV^. These discussions contain the zenith angle depen- 
dence's of sub-GeV and/or multi-GeV Kamiokande data. The zenith angle dependence's in 
recent SuperKamiokande atmospheric experiments 0,0] have decreasing character not only 
for multi-GeV experiment but also for sub-GeV one. We have analyzed these dependence's 
partially |]12| and obtained a result that large angle solution for sin^ 2^12 is favorable. In this 
paper, we will analyze these dependence's thoroughly. 



II. NEUTRINO OSCILLATION IN THREE-FLAVOR NEUTRINOS 



Weak currents for the interactions producing and absorbing neutrinos are described as 



o,/3=l 



(1) 



where la {h = e, I2 = /i, h = t) represents the lepton flavor, the neutrino mass eigenstate 
and U is the lepton mixing matrix. U is the unitary matrix corresponding to the CKM matrix 
^CKM quarks defined by 



U = UiUl 



(2) 



where the unitary matrices Ui and 11^ transform mass matrices M' for charged leptons and 
M'^ for neutrinos to diagonal mass matrices as 



UiM^Uf'^ = diag [me, m^,mr], 
UyM^U'^ = diag[mi,m2,m3]. 
We present the unitary matrix neglecting the CP violation phases as 



(3) 



JJ _ g«6'23A7gi6l3A5gi6'l2A2 



V 



'12'-^23 

■' a'^ 
'12'' 23 



12'-13 ^12'-13 ■'13 

_ „u „u „v „v „v „u „u „u „v 

12''23''13 ''12''23 ■'12''23''13 ''23''13 



-^12'-23''13 



COS^,^ 



where Aj's are Gell-Mann matrices. 
The probabilities for transitions vi^ 



"'-12"' 23 

sin Q. 



'12'-23''13 '-23'-13 



(4) 



vi are written as 



S23 + p. 



31 



s- 



p23 



-25, 



"'Q "'^ 

w„a - 2Ul,) + 2{-Ul,Ul, + UlM,2 + Ul^Uf.)^ 



31, 



-25,„,,(1 - 2UfJ + 2{Ul,Ul, + Ul^Uf. - Ul^, 



,31 



-25^.^,(1 



+ 2{Ul,U'-Ul,U?, + Ul,U? 



(5) 
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where Sij is the term representing the neutrino oscillation; 

Sij = sm 1.27— ^L, 



(6) 



in which Am?- 

''J 



m7 



'J I' 



10 ^°eV^ is obtained |T9[, and for the atmospheric experiments, in 



E and L are measured in units eV^, GeV and km, respectively. 
The values of neutrino masses are not known precisely, but two mass parameters are 
necessary to account for the solar neutrino experiments, in which mass parameter Am^ is 
10"^ - lO^^eV^ or ~ 
which Am^ is 10^^ 

is that lower two neutrino masses in three neutrinos are very close and another one is rather 
far away from them. Then we assume that three neutrino masses have such a mass hierarchy 
as 



lO^W^ 19 . Then it seems that the most economical mass hierarchy 



(7) 



In this mass hierarchy Eq. (|^, Am^2 ^ ^'^Is — ^"^13 ' expression Eq. (5) for the 



transition probabilities P{jyi^ 

1 - 2(1 



Ui ) are rewritten as 



1 - 2(1 - 2\Jl, - Ul, - t/i + 



t2 



P{iye 
P{ye 
P{y.. 



2(1 - 2f/4 - Ut, - + Ut,)S,2 - 4f/4(l - f/4)S23, 
2{Ul,Ul + U%Ul, - U%Ul,)S,, + AUl,U%t 



^J■J 



23, 



(8) 



III. NUMERICAL ANALYSES OF NEUTRINO OSCILLATIONS 



A. Atmospheric neutrinos 

The evidence for an anomaly in atmospheric neutrino experiments was pointed out by 



Kamiokande Collaboration and 1MB Collaboration |^ using the water-Cherencov 
experiments. More recently, SuperKamiokande Collaboration reports the more 

precise results on anomaly in atmospheric neutrino. The double ratios i?(/i/e) = 
-Rexpt(/^/e)/-RMc(/^/e) obtained are summarized as 

, 0.60+°-°^ ±0.05 for Kamiokande (sub-GeV) [3], 
i?(/i/e) = { ■ (9a) 

0.57+°°^ ±0.07 for Kamiokande (multi-GeV) H, 
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R[li/e) = 0.54 ± 0.02 ± 0.07 for 1MB |2|, (9b) 



0.63 ±0.03 ±0.05 for SuperKamiokande (sub-GeV) 0, 
i?(/i/e) = { (9c) 

0.60 ± 0.06 ± 0.07 for SuperKamiokande (multi-GeV) §. 

In the above data, sub-GeV experiments detect the visible-energy less than 1.33GeV. The 
zenith angle dependence's of double ratio i?(/i/e) obtained in SuperKamiokande are shown 
in Fig. 1 [0. The sub-GeV data has decreasing character to increasing of zenith angle 
in contrast to the zenith angle dependence of Kamiokande data which has no decreasing 
character [^. The multi-GeV data has also decreasing character to increasing of zenith 
angle similar to that of Kamiokande [^] . 

The ratios i?expt(Ai/e) and i?Mc(/^/e) are defined as 

/texptl/i/ej j2aIeeiEe)ae{E,,E,)F^{E,,9)Piu^ ^ uMEedE^de ' ^ 
„ . Je,{E^)a,iE,,E^)F^{E,,9)dE,dE,de 

nucWe) J ,^(^E,)^^{E,^E,)F,{E,,e)dEJE,de ' ^'^^^ 

where the summation J2a ci^re taken in /z, e neutrino, and in these expressions antineutrino 
processes are contained. ep^Ef^) is the detection efficiency of the detector for /3-type charged 
lepton with energy Ej^, is the differential cross section of and Fa{Eu, 6) is the incident 
Va flux with energy E^ and zenith angle 9. P^Ua i^p) is the transition probability Eq. (8) 
and it depends on the en ergy E^, and the distance L which depends on zenith angle 6 as 
L = \J {r + hY — sin^ 9 — rcos6, where r is the radius of the Earth and h is the altitude 
of production point of atmospheric neutrino. The zenith angle dependence's of double ratio 
are defined as the equation similar to the Eq. (10) except for the no integration on zenith 
angle 9. 

Although information of Fa{E^,9) etc. are given in Refs. p5H27||, we use the MC pre- 
dictions for fa{Ey,9) = /e/,(e)(-E^{e))a^{e)(-Ei.,£'^(e))F«(£'^,6')(iE^(e) in Ref. [|] for sub-GeV 
experiment and Ref. for multi-GeV experiment. Explicit E^ dependence of fa{Ea,9a) 
are shown in Appendix B of Ref. [jl2[. Since Pii^a ~^ ^pi) and P{va ^e) are the functions 
of (Am^2) ^^23' ^12) ^13) ^23; -^5 -S'); the double ratio R{fi/e) which is integrated in neutrino 
energy E and zenith angle 9 (related to L) is a function of (Am^g? "^^23' ^i2) ^is? ^23)- The 
zenith angle dependence of double ratios is a function of (Am^2) ■^^23) ^12 > ^13 > ^23, 9). 



In previous paper [0, we estimated the double ratio R{fj,/e) fixing the parameters 



(Am^2 5 sin^ 2^) which have been determined in solar neutrino experimental data [pl|-|2^ 



using the three-fiavor neutrino analyses |I^p!2| , as Amfg = 3 x 10~^eV^, sin^ 26'i2 = 0.7 
which corresponds to large angle solution and Am^2 — 10~^eV^, sin^ 26'i2 = 0.005 which 
corresponds to small angle solution. We showed these plots in Fig. 2. Contour lines 
correspond to the upper and lower values of R{fi/e) in SuperKamiokande data. We showed 
the plots of sub-GeV experiment in Figs. 2(a)-(d) and plots of multi-GeV one in Figs. 2(e)- 
(f), and in these figures solid lines denote the large angle solution plots and dotted lines 
the small angle solution plots. In Figs. 2(e)-(h), dotted lines are close to solid lines. Same 
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analysis has been performed by the third paper in Ref. jl^ and our results obtained are 
similar to those. 

In Fig. 1^, we showed the contour plots of taken from the zenith angle dependence 
data (values drawn from Fig. 1) for double ratio R{fj,/e) on tan^ ^13 — tan^6'23 plane for 
various Aml^. Contour curves denoted by dotted, thick solid and thin solid hne correspond 
to the values of = 7.78 (90% C.L. for 4DF), 9.49 (95% C.L. for 4DF) and 13.3 (99% 
C.L. for 4DF), respectively. We showed the plots of sub-GeV experiment in Fig. |^(a)- 
(d) and plots of multi-GeV one in Fig. ^(e)-(h). These plots are taken for large angle 
solutions of sin^ 2612 because for small angle solution cannot be less than 12 in sub- 
GeV experiments. This is seen by the decreasing character of zenith angle dependence 
in sub-GeV data. We showed the predicted zenith angle dependence curves of R{fi/e) on 
the graph of zenith angle dependence data in Fig. 1. Fig. 1(a) represents the sub-GeV 
experimental case on parameters {Am^^ = O.leV^, 613 = 2.4°, 623 = 25°). In this figure, 
solid curve denotes the large angle solution case and dotted curve the small angle solution 
case. From this figure, we can understand the small angle solution case do not give the 
small x^- It can be said that the decreasing character of sub-GeV data is caused by the 
large fi — e mixing. In Fig. 1(b), we showed the multi-GeV experimental case on parameters 
(Amig = 0.1eV^ ^13 = 4°, 623 = 45°). In this case, there is no difference between large 
angle solution and small one. Decreasing character in multi-GeV data is explained by the 
large /i — r mixing. 

In Fig. 4, we showed the contour plots of of sub-GeV plus multi-GeV experiments. 
These are taken for large angle solution. Contour curves denoted by dotted , thick solid 
and thin solid line correspond to the values of x^ 14.7(90% C.L. for 9DF), 16.9(95% C.L. 
for 9DF) and 21.7(99% C.L. for 9DF), respectively. Thus, if we claim that both sub- 
GeV data and multi-GeV one should be explained simultaneously, the allowed region on 
parameters become narrow. The best fit by for sub-GeV plus multi-GeV data is obtained 
on parameters (Am23 = O.leV^, 6*13 = 33°, 623 = 45°) on which x^ is 13.1. Although 
the allowed region on parameters is obtained from the atmospheric neutrino experiments, 
all regions are not accepted because the terrestrial neutrino experiments seem unfavorable 
for the neutrino oscillation. In next subsection, we will discuss the terrestrial neutrino 
experiments, thoroughly. 



B. Terrestrial neutrinos 



Although we have study the terrestrial neutrino experiments in previous paper []T2 
we did not include the Chooz experiment . 
the short baseline experiments, there are E531 
experiments searching for Ur appearance in u^. 
NOMAD experiments; 



We include it in this present analysis. In 
§, CHORUS and NOMAD § accelerator 
We used the data of E531, CHORUS and 



P(z^M ^ z^r) < 2 X 10"^ (90% C.L.), 



(11) 
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L/E ~ 0.02. 



For the experiments searching for z/^ and z/^ — > Ve oscillations, we used the experiments 
E776 [0, KARMEN § and LSND § accelerator experiments; 

P(z/^ ^ z/e) < 3 X 10"^ (90% C.L.), E776 (12a) 

L = 1km, E ~ IGeV, 
^(z^M ^ '^e) < 3.1 X 10-3 (90% C.L.), KARMEN (12b) 

L = 17.5m, E < 50MeV, 
i^(z>M ^ '^e) = 3.4+2-0 ± 0.7 X 10-^ LSND (12c) 

L = 30m, E ~ 36 - 60MeV. 



We analyzed the Bugey experiment using nuclear power reactor |10| searching for the dis- 



appearance of z/e. In this paper, we add the Chooz experiment |TT| . 

1 - P(z/e ^ z/e) < 10^2 (90% C.L.), Bugey (13a) 

L = 15, 40, 95m, P ~ 1 - 6MeV. 
1 - P(z/e ^ z/e) < 0.9 X 10"^ (90% C.L.), Chooz (13b) 

L = 1km, E ~ 3MeV. 

We show the contour plots of the allowed regions on (tan^ ^13, tan^ ^23) plane determined 
by the probability P expressed in Eq. (^ and above experimental data Eqs. (pAf), (|I2[) 
(except LSND data ( |12c| )) and (p!3D in Fig. ||. Allowed regions are left and right hand sides 
surrounded by curves. Curves represent the boundary of 90 % C.L. of P. We fixed the 
parameters Am^2 and 612 as Am^2 = lO-^eV^ and sin^ 26'i2 = 0.8, and the parameter Aml^ 
to be various values from O.OOleV^ to leV^. Although we fix the parameters Am^2 and 
612 as Am^2 = lO-^eV^ and sin^ 2^12 = 0.01, the results are not changed because the 
probability P of terrestrial neutrino is insensitive to the e — /i mixing parameters Amf 2 and 
sin2 26'i2. Dotted lines show the allowed regions restricted by the LSND data Eq. ( |12c[ ). 
From these results, we see that the allowed regions on {613, 6*23) is very restricted by LSND 
data. Comparing the result obtained previous analysis with this, there is observed that 
on mass parameter Am|3 ~ 10"^ — 2 x 10"^, allowed region for tan^ ^13 in present result is 
rather restrictive than the one in previous result. 



C. Allowed regions of mixing parameters 

We here discuss the allowed regions for mixing parameters (Amf2, sin^ 26'i2, Am23, ^13, 
^23) satisfying atmospheric neutrino and terrestrial neutrino experiments. Am^2 and sin^ 2^12 
have been determined in analyses considering MSW effects of the solar neutrino experiments 
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(Ami2,sm^ 26*12) 

(4 X 10-6 - 7 X 10-^eV^ 0.6 - 0.9), large angle solution 

(3 X 10-6 - 1.2 X 10-^eV^ 0.003 - 0.01). small angle solution 
for = 0° - 20° 

Observing the allowed regions obtained in terrestrial neutrino experiments (Fig. P) and 
the allowed regions obtained in atmospheric neutrino experiments including sub-GeV and 
multi-GeV data but not including zenith angle dependence's (Fig. ^), we obtain the allowed 
regions satisfying these experiments. 

large angle solution small angle solution 

for Am23 = lOeV^ no allowed region no allowed region (15a) 

for Am^g = leV^ (^13 < 4°, 24° < ^23 < 62°) (^13 < 4°, 24° < ^23 < 66°) (15b) 

for Am^g = O.leV^ (^13 < 3°, 27° < ^23 < 63°) (^13 < 3°, 27° < ^23 < 63°) (15c) 

for Am^g = O.OleV^ (^13 < 13°, ^23 ~ 45°) (^13 < 13°, ^23 ~ 45°) (15d) 

for Am23 ~ O.OOleV^ no allowed region no allowed region (15e) 

We next show the allowed regions satisfying the terrestrial experiment data (Fig. 5) and 
zenith angle dependence's of atmospheric experiments (Fig. 3). We showed the allowed 
regions for sub-GeV data and multi-GeV data, separately. 

For sub-GeV case: allowed region are obtained for large angle solution. 

for Am23 = lOeV^ no allowed region (16a) 

for Am^g = leV^ (^13 < 4°, 18° < ^23 < 32°) (16b) 

for Am^g = O.leV^ (^13 < 3°, 18° < ^23 < 33°) (16c) 

for Am^g = O.OleV^ (^13 < 13°, 17° < ^23 < 32°) (16d) 

for Am^g = O.OOleV^ (^13 < 23°, 21° < ^23 < 35°) (16e) 

For multi-GeV case: allowed region are obtained for large and small angle solutions. 

for Am23 = lOeV^ no allowed region (17a) 

for Am^g = leV^ (^13 < 4°, 27° < ^23 < 63°) (17b) 

for Am^g = O.leV^ (^13 < 3°, 28° < ^23 < 61°) (17c) 

for Am^g = O.OleV^ (^13 < 13°, 35° < ^23 < 55°) (17d) 

for Am^g = O.OOleV^ (11° < ^13 < 23°, 45° < ^23 < 61°) (17e) 

We finally show the allowed regions satisfying the terrestrial experimental data (Fig. 5) 
and atmospheric experiments including the zenith angle dependence's for both sub-GeV and 
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multi-GeV data (Fig. 4). 



For sub-GeV plus multi-GeV case: allowed region are obtained for large angle solution 
and Aml^ = 0.08 - 2eVl 

for Am^g = 2eV^ (^13 < 4°, 623 ~ 30°) (18a) 

for Am^g = leV^ (^13 < 4°, 27° < 623 < 32°) (18b) 

for Aml^ = O.leV^ (^13 < 3°, 28° < 623 < 33°) (18c) 

for Amis = 0.08eV2 (^13 < 3°, 623 ~ 30°) (18d) 

These regions are very restricted than atmospheric experiments regions Eqs. (15) obtained 
by not considering the zenith angle dependence's. Although the LSND data is very restrictive 
experiment, there exists allowed region satisfying this experiment at the values of parameters 
(Am^g ~ leV^ = 3 - 4°, 623 = 27- 32°) with large angle solution. 



IV. DISCUSSIONS 



We analyzed atmospheric neutrino experiments of SuperKamiokande [JI],0 using the 
three-flavor neutrino framework with the mass hierarchy mi ^ m2 -C ^3 and obtained 



n2, 



the allowed regions of the parameters (Am^g, sin^ 2^ 
the terrestrial neutrino experiments including Chooz data |TT 



Amis, ^13 



6^23). We also analyzed 
From these neutrino ex- 
periments (not considering the zenith angle dependence) and the results obtained in solar 
neutrino experiments we obtained the rather wide allowed regions denoted in Eqs. 
(15). When zenith angle dependence's are considered, the allowed regions are restricted as 



follows: (Am 
and (6*13 < 4' 



2 

12 



= 4 X IQ-^ - 7 X 10 
27° < 023 < 32 



sm 



29 



12 



for Am23 



leV 



0.6 



'13 



- 0.9, 

< 3°, 



Amis ■ 
28° < 



= 0.08 

^23 < 



- 2eV^), 
33°) for 



Amis 



O.leV . If we include LSND experiment in terrestrial experiments, allowed regions 



are restricted to (Am^g ~ leV^ ^13 = 3 - 4°, ^23 = 27 - 32° 



Finally, we present the neutrino mixing matrix Eq. 
solutions Eqs. (18); 



(2) numerically for the allowed 



^ 0.81^ 0.90 0.43 0.58 0.0 ^0.07^ 



U 



-0.36 



-0.55 0.66 ^ 0.80 0.45 ^ 0.54 



0.14 ^ 0.32 



-0.37 ^ -0.51 0.84 ^ 0.89 



(19) 
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FIGURES 

FIG. 1. The zenith angle dependence of R{n/e) in atmospheric neutrino experiment. Ex- 
perimental data is referred to SupcrKamiokande (Ref. [2]). Fig. (a) represents the sub-GeV 
experimental case: solid curve is calculated on typical parameters for large angle solution (Amf2 = 
3 X lO-^eV^, sin2 26'i2 = 0.7, Am23 = 0-2eV'^, = 4°, 023 = 25°) and dotted curve on small 
angle solution (Aml^ = lO^^eV^ sin^ 26*12 = 0.007, Am23 = 0.2eV^ 6*13 = 4°, 6*23 = 25°). 
Fig. (b) represents the multi-GeV experimental case: solid curve is calculated on parameters 
(^13 = 4°, 023 = 30°) for large angle solution and dotted curve on parameters (^13 = 4°, ^23 = 30°) 
for small angle solution. 

FIG. 2. The plots of allowed regions determined by atmospheric neutrino data Eq. (9c) of 
R{lJ./e) for various values of Am^s on tan^ ^13 — tan^ ^23 -Figs, (a)-(d) show the plots of sub-GeV 
experiments and Figs, (e)-(h) the ones of multi-GeV experiments. We fix the parameters (Amf2, 
sin^20) to be (3 x 10~^eV^, 0.7) corresponding to the large angle solution (solid lines) and 
(10~^eV^, 0.005) corresponding to the small angle solution (dotted lines). 

FIG. 3. The plots of allowed regions on tan^ ^13 — tan^ ^23 plane determined by the zenith 
angle dependence's of SuperKamiokande data using the analysis. Figs, (a)-(d) shows the plots 
of sub-GeV experiments for large angle solution of sin^ 2^12 and Figs, (e)-(h) the ones of multi-GeV 
experiments for large and small angle solutions. Solid thin, solid thick and dotted curves denote 
the regions allowed in 99%, 95% and 90% C.L., respectively. 

FIG. 4. The plots of allowed regions on tan^ ^13 — tan^ ^23 plane determined by the sub-GeV 
and multi-GeV zenith angle dependence's of SuperKamiokande data using the analysis. These 
are taken for large angle solution. Solid thin, solid thick and dotted curves denote the regions 
allowed in 99%, 95% and 90% C.L., respectively. 
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FIG. 5. The plots of allowed regions on tan^ ^13 — tan^ 823 plane determined by P of terrestrial 
z^/i ^T, ^fj, z^ej i^fj, and Vf, experiments. Allowed regions are left and right hand 

sides surround by curves and/or lines. Curves and lines represent the boundary of 90 % C.L. of P. 
Amf2 and sin^ 2^12 are fixed as 10~^eV^ and 0.8, respectively. Am23 fixed to leV^(Fig. 5(a)), 
0.1eV2(Fig. 5(b)), 0.01eV2(Fig. 5(c)) and 0.001eV2(Fig. 5(d)). Dotted lines show the allowed 
regions determined by LSND data. 



13 



1.2- 





o 



DC 0.8- 



°1 



Super-Kamiokande Preliminary 
325.8 days Sub-GeV 1 ring 




large angle sol. 
small angle sol. 





COS0 



Fig. 1(a) 



1.2- 



Super-Kamiokande Preliminary 
325.8 days Multi-GeV FC+PC 





o 



DC 0.8- 




< 

^2 0.4 



°-1 



large angle sol. 
small angle sol, 





COS0 



Fig. 1(b) 




Fig. 2 




Fig. 2 




90%(CL) 
95%(CL) (4DF) 
99%(CL) 




90%(CL) 
95%(CL) (4DF) 
99%(CL) ^'"23^ 



AmoVO.leV^ 



tan 01 



tan 9,. 



(a) 



(b) 




90%(CL) 
95%(CL) (4DF) 
99%(CL) 



<^ o 




90%(CL) 
95%(CL) (4DF) 
99%(CL) "'"23 



Am2^-.=0.001eV^ 



tan 01 



tan 9i- 



(C) 



(d) 



Fig. 3 




90%(CL) 
95%(CL) (4DF) 
99%(CL) 




90%(CL) 
95%(CL) (4DF) 
99%(CL) 



tan 01 



tan Qr 



(e) 



(f) 



90%(CL) 
95%(CL) (4DF) 
99%(CL) 



<^ o 




90%(CL) 
95%(CL) (4DF) 
99%(CL) "'"23 



Am2^-.=0.001eV^ 



tan 01 



tan 9i- 



(g) 



(h) 



Fig. 3 



90%(CL) 
95%(CL) (9DF) 
99%(CL) 



tan 01 



90%(CL) 
95%(CL) (9DF) 
99%(CL) 



AmoVO.leV^ 



10 10 10 10 10 10 10 



tan ^9] 



'13 



(a) 



(b) 



CP 



90%(CL) 
95%(CL) (9DF) 
99%(CL) 



Am|,=0.01eV^ 



tan 01 



<^ o 



90%(CL) 
95%(CL) (9DF) 
99%(CL) 



Am2^-.=0.001eV^ 



10 10 10 " 10 " ^ 10 10 10 



tan ^9] 



'13 



(c) 



(d) 



Fig. 4 





tan 01 



tan Qi 



(a) 



(b) 





tan 01 



tan 9i- 



(C) 



(d) 



Fig. 5 



